The nature of slow dynamics in a minimal model of frustration-limited domains 
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We present simulation results for the dynamics of a schematic model based on the frustration- 
limited domain picture of glass-forming liquids. These results are compared with approximate 
theoretical predictions analogous to those commonly used for supercooled liquid dynamics. Although 
model relaxation times increase by several orders of magnitude in a non-Arrhenius manner as a 
microphase separation transition is approached, the slow relaxation is in many ways dissimilar to 
that of a liquid. In particular, structural relaxation is nearly exponential in time at each wave 
vector, indicating that the mode coupling effects dominating liquid relaxation are comparatively 
weak within this model. Relaxation properties of the model are instead well reproduced by the 
simplest dynamical extension of a static Hartree approximation. This approach is qualitatively 
accurate even for temperatures at which the mode coupling approximation predicts loss of ergodicity. 
These results suggest that the thermodynamically disordered phase of such a minimal model poorly 
caricatures the slow dynamics of a liquid near its glass transition. 
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I. INTRODUCTION 

Several microscopic scenarios have been proposed as 
underlying mechanisms for dynamical arrest in super- 
cooled liquids 0, S H S 0. Because this dramatic 
slowing down is accompanied by the onset of dynam- 
ical heterogeneity [(J 0, 0, a promising candidate model 
should account for the spontaneous segregation of rapidly 
relaxing and slowly relaxing domains. The theory of 
frustration-limited domains has been developed with this 
condition in mindj^, The basic units of this the- 

ory are microscopic regions of low internal energy whose 
spatial extent is limited by long-ranged interactions or 
constraints. Nelson and coworkers suggested that such 
domains form in simple "atomic" liquids (such as metal- 
lic glass- formers') . For small clusters of particles, 
icosahedral arrangements are energetically preferred over 
closest-packed configurations representative of the crys- 
talline state^^. It is thus argued that a supercooled 
atomic liquid is rich in low-energy icosahedral clusters 
which, for geometric reasons, cannot extend indefinitely. 
In this case, frustration is a consequence of the vanishing 
curvature of Euclidean space. One can imagine that the 
nature of local order and source of frustration are some- 
what different for other, more complicated materials. 

The frustration-limited domain theory has several at- 
tractive features. By associating relaxation kinetics with 
the interfacial area of domain walls, it predicts a crossover 
in the temperature dependence of structural rearrange- 
ment times from Arrhenius to a super- Arrhcnius form. 
Such a crossover is prominent in experiments with frag- 
ile glass-formers |l3| , Further, the theoretical scaling ex- 
ponent for the asymptotic super- Arrhenius temperature 
dependence appears to be consistent with experimen- 
tal data. Finally, the notion of frustration-limited do- 
mains resonates with the picture of heterogeneous dy- 



namics that has emerged from experiments 0, 0, @ and 
simulations [H [H EH • 

While local structure within low- 
energy clusters may be effectively frozen on a molecular 
time scale, relaxation can be facile at the strained in- 
terfaces between domains. This argument is the essence 
of Stillinger's "tear and repair" picture of shear flow in 
fragile liquids 0. 

Although these dynamical predictions of frustration- 
limited domain theory are suggestive, they are funda- 
mentally thermodynamic in nature and thus indirect. In 
order to make such arguments precise, Grousson et al. 
have recently focused on the explicit dynamics of min- 
imal models exhibiting frustration- limited domains |l8j. 
Specifically, they have simulated stochastic dynamics of 
several classical spin models that pit short-ranged, fer- 
romagnetic interactions against long-ranged, antiferro- 
magnetic interactions. In addition to confirming super- 
Arrhenius relaxation for these models at low tempera- 
tures, they have demonstrated that the "fragility" of the 
dynamics (i.e., the degree of deviation from Arrhenius 
form) varies continuously with the relative strength of 
the long-ranged frustration. Quite recently, Grousson et 
al. have performed dynamical mode coupling calcula- 
tions for the same models, but no direct comparisons to 
their earlier simulations were made[T^. In this paper, we 
perform similar calculations that are compared directly 
to numerical simulations. 

Schmalian and coworkers have also argued that such a 
Coulomb-frustrated ferromagnet should display essential 
features of glassy dynamics pol EH E5 | . In their analysis, 
it is a proliferation of metastable states that drives vit- 
rification. With the aid of replica mean field theory, this 
perspective predicts a scaling of fragility with frustra- 
tion str eng th that agrees well with the simulation results 
of Ref. 18] . But like the theory of frustration-limited 
domains, this analysis is thermodynamic in nature, re- 
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lying on an assumed correspondence between particular 
subensembles of high free energy and genuine dynamical 
bottlenecks. 

This paper addresses the extent to which the slow re- 
laxation in such simplified models truly resembles that of 
molecular glass-forming liquids. For this purpose we in- 
vestigate in detail the dynamics of a model closely related 
to those studied in Refs.fl8| and j2jj. The model and dy- 
namical propagation rules we consider, which are free of 
artificially quenched disorder and kinetic constraints, are 
described in Sec. 2. For several values of the frustration 
strength, we compare the time dependence of spin corre- 
lations computed in simulations with those predicted by 
approximate theoretical approaches. 

Two self-consistent dynamical equations obtained from 
theory are discussed in Sec. 3. They correspond to re- 
summations of an exact, infinite diagrammatic series for 
time correlations. The first resummation, yielding expo- 
nential relaxation at each wave vector, is a direct dynam- 
ical generalization of Brazovskii's static result for this 
class of models 23] . Calculations based on this straight- 
forward approach agree remarkably well with simulation 
results, described in Sec. 4, even for temperatures ap- 
proaching a thermodynamic transition to a fully ordered 
state. The second resummation is formally analogous to 
the idealized mode coupling theory of liquids. As such, 
it predicts loss of ergodicity at finite temperature. By 
contrast, we find no evidence of nonergodic behavior or 
two-step relaxation for simulated disordered states of this 
model. 

Surprisingly, then, the Hartree approach is the more 
accurate approximation for slowly relaxing disordered 
states of the model system. Detailed comparison of simu- 
lation and theory confirms that sluggishness indeed arises 
from de Gennes narrowing (i.e., from dramatic changes 
in static correlations), in contrast to complex dynamical 
mechanisms such as mode coupling. In this respect, the 
minimal model we have studied does not capture impor- 
tant aspects of supercooled liquid dynamics, specifically 
the intermediate time plateau and long time stretching 
of dynamical correlators. Implications of this result are 
discussed in Sec. 5, along with issues related to fragility 
and local conservation of magnetization. In section 6 we 
conclude. 



II. MODEL 

We consider fluctuations of a field 4>(r) at position r in 
three dimensions, with energy |23L I24I Eq 



W(r)] 



dr 



<Mr) (r + fc - 2 (V 2 + fc 2 ) 2 ) 



(r) + ^V) 



(1) 



Here, the energy scale characterizes typical fluctua- 
tions of a surrounding heat bath. The physical meaning 
of the field 4> may be somewhat abstract in the context 



of supercooled liquids, for example representing the de- 
gree of a particular local packing symmetry. The ap- 
plication of Eq- CD to diblock copolymer melts is more 
intuitive[2H l28l l2^ | . In this case, <fi represents the local 
excess number density of one monomer type, and r de- 
scribes the preferential affinity of monomers for others of 
the same type. Because our interest in this model is moti- 
vated by the work in Ref.|l8||. we will imagine that <fi sim- 
ply represents a coarse-grained, scalar spin density. Here, 
t is a dimcnsionless temperature measuring the distance 
from an underlying critical temperature when fco = 0. 
For any physical interpretation, the wave vector ^ / 
characterizes long-ranged order of a low-temperature, mi- 
crophase separated state. The coefficient A multiplying 
the c/> 4 term in Eq. ^will later be used to order terms in 
perturbation series, although in calculations its numeri- 
cal value will be of order unity. 

For wave vectors near fco, the action in Eq. ^ cor ~ 
responds to that studied by Schmalian and Wolyncs 
20, 2lll22| and (in a hard-spin lattice version) by Grous- 
son et q?.|18j. In their work, frustration is explicit in 
competing interactions of square-gradient (J|V0(r)| 2 ) 
and Coulomb (Q 4>(r)(j>(r')/\r — r'|) forms. The rela- 
tive strengths of these interactions determine the peri- 
odicity of the ground state, fco oc (Q/J) 1 ^ 4 , in which 
spin-up and spin-down domains alternate in stripes or 
lamellae 26| . In the model defined by Eq. ^ this frustra- 
tion is instead implicit in nonzero fco, but has the same 
physical effect. Namely, homogeneous domains are ener- 
getically favored at small length scales, while net mag- 
netization is effectively constrained to vanish at larger 
length scales. In the context of diblock copolymers, this 
effective constraint reflects the stoichiometry imposed by 
polymer connectivity [271 l28t |29| . 

In two dimensions and higher, the presence of nonzero 
fco in Eq. ^ff as a subtle but profound effect on the ther- 
modynamics of the paramagnetic state. Specifically, the 
large entropy of fluctuations near |k| = fco significantly 
reduces the free energy of the disordered phase. Within 
a Hartree approximation, this contribution is sufficient 
to make the paramagnetic state stable or metastable for 
all finite r. As a consequence, the transition to a phase 
with long-ranged order is first (rather than second) order 
and occurs at a temperature r tr < 0. This effect was 
first recognized by Brazovskii|23j and has been summa- 
rized lucidly by Binder and Fredrickson|29j . Its quali- 
tative features have been subseque ntly confirmed in ex- 
periments with diblock copolymers[30j. Because, in this 
picture, statistics of the disordered state are dominated 
by fluctuations near fco, we expect the model of Eq. ^to 
belong to the same universality class as those of Refs.|l8j| 
and pbl I2H |2^. Later, we will demonstrate that slow 
dynamics of this state are dominated by the very same 
fluctuations. 

The action in Eq. ^ is n °t a true Hamiltonian, and 
thus has no intrinsic dynamics. We consider two com- 
monly used stochastic propagation rules which generate a 
canonical ensemble of fluctuations consistent with Eq. ^ 



3 



Attention will be primarily focused on a simple Langevin 
equation, 



at 



6</>(r) 



+ ??(r,i), 



(2) 



where rj(r, t) is a random force whose statistics are Gaus- 
sian, and 



(7?(r, t)f](r',t')) = 2p- 1 S{r - r')S(t - t'). 



(3) 



In Eq. |3 angled brackets denote an average over all 
possible realizations of the random force. The above 
equation of motion, along with the energetics of Eq. ^ 
and the statistics of Eq. |3| has been studied previously, 
most notably in the context of nucleation and nonequi- 
librium pattern formation following a rapid quench to 

our knowledge, the detailed 
equilibrium dynamics of the paramagnetic phase very 
close to the transition (i.e., r > r tr ) have not until now 
been fully explored. 

The dynamics generated by Eq. |2 do not conserve the 
field <Xr, £)[33|]. Because the slow relaxation of super- 
cooled liquids results in part from the conservation of 
hydrodynamic densities, this feature of Eq. [21 may be 
somewhat troubling (particularly in the context of di- 
block copolymers, in which the number density is clearly 
conserved). For this reason, we consider a second form 
of dynamics that conserves </>(r, t) by construction. Tra- 
jectories of this dynamics are chains of microstates gen- 
erated by a Metropolis Monte Carlo algorithm. In detail, 
a random displacement of the field </>(r) is attempted at 
discrete time steps, and is accepted with probability 



P acc = mm[l,exp(-/3AH)] 



(4) 



where A7i is the change in energy produced by the dis- 
placement. Local conservation of the field is achieved by 
restricting the choice of random displacements to those 
which do not alter the local net magnetization. Further 
details of this procedure will be described in Sec. 4. Sim- 
ulation results presented in that section demonstrate that 
the decay of spin correlations produced by Eqs.[21and0] 
are nearly identical, within an arbitrary rescaling of time 
in the Monte Carlo chain of states. The physical meaning 
of this fact will be discussed in Sec. 5. 



III. THEORY 

In this section we discuss two approximations for re- 
laxation of the field (/>(r,i). Specifically, we derive closed 
equations of motion for the correlation function 



Ck(t-i') = <0 k (i)0-k(i')), 



(5) 



where Fourier components of the field are defined in the 
standard way: 



Mt) = J dr0(r,t)e lk r . 




(6) 



FIG. 1: Diagrams representing low-order terms in a pertur- 
bation series for the dynamics generated by Eq. 2. Lines rep- 
resent instantaneous correlations of the field at two points in 
space, while vertices represent interactions. 



The first equation of motion is linear in (f>k(t) and resem- 
bles phenomenological theories of high temperature liq- 
uid state dynamics, such as the wave vector dependent 
viscoelastic theory |34j. The second is nonlinear and has 
the form of the idealized mode-coupling approximation 
to the dynamics of density fluctuations in supercooled 
liquids 0. This equation contains the feedback mecha- 
nism responsible for the interruption of particle diffusion 
at intermediate time scales due to constraints imposed 
by slowly reorganizing local environments (the "cage" ef- 
fect). 

Treating A as a perturbation parameter, the solution to 
Eq. |2 may be written as an infinite series of terms, each 
representing a collection of field interactions and periods 
of free propagation. As a result, Ck(t) and its associ- 
ated response function, Gk(t) = —(3 dCk(t)/dt can be ex- 
panded in powers of A. (We focus exclusively on the por- 
tion of the phase diagram in which dynamics are ergodic, 
so that Ck(£) and Gk(i) are related by the fluctuation- 
dissipation theorem.) Fig. ^ shows diagrammatic repre- 
sentations of the first two terms in the series for C\(t). 
The development of this expansion, as well as the par- 
tial series resummations underlying our approximations, 
have been discussed thoroughly in the context of other 
models. We will describe physically significant highlights 
of the procedure and refer the reader to Refs.[l!j[35| for 
details. 

A linear equation of motion for Gk(i) results from sum- 
ming only terms in the series whose diagrams have the 
basic topology shown in Fig.^a). In the irreducible seg- 
ments of these "tadpole" diagrams, all interactions coin- 
cide in time. Consequently, such a summation renormal- 
izes only the static portion of the basic tadpole diagram 
of Fig.^a). The dynamics predicted by this Hartree re- 
summation scheme are identically those of a variationally 
optimized harmonic reference system [3^. They are thus 
obtained more directly by assuming Gaussian statistics 
for 4>k(t). Specifically, we multiply the Fourier transform 
of Eq. |3 by 0_k(O) and average over the noise history, 



4 



yielding 
dt 



+ k Q 2 (k 2 -k 2 ) 2 }Cu(t) + 



A 



- ]T (^k/(t)^(t)^k_k'-k«(*)0-k(O)>. (7) 



k'.k" 



Eq. [7| is the first member of a complicated hierarchy of 
equations relating multi-point fluctuations to correlations 
of higher order. But if we assume that ^(t) is a Gaussian 
random variable, the hierarchy closes immediately: 



dt 



(8) 



Here, the renormalized mass th that appears in the 
expression for the structure factor /ij? = 1/Cjr(0) = 
Th + k^ 2 (k 2 — k 2 ,) 2 is determined self-consistently by 



E^(o). 



(9) 



This result is precisely Brazovskii's static 
approximation!^. The relaxation described by 
Eq. [SI while simply exponential, occurs with rates 
that are significantly renormalized by the entropy of 
fluctuations near |k| = fco. 

More elaborate, nonlinear approximations for Ck(i) re- 
sult from incorporating diagrams with more complicated 
topologies |l9t l35| . The mode coupling approximation 
(MCA) is an example, including diagrams with the "sun- 
set" shape of Fig. H(b) . Summing all terms that renor- 
malize the propagators (but not the vertices) of the basic 
sunset diagram yields the MCA. Because these contribu- 
tions involve more than one unique time variable, they 
are capable in principle of capturing nontrivial memory 
effects. As shown in Ref.[35|, the self-consistent result of 
this resummation is 



dc MCA (t) t 



= -Mk* CA Cf CA W + 



dt 
X 2 p '* 
6 



/ dt V [c™ CA (t-t')c™F A (t-t') 
Ja k/,k« 

xc^iu»(t-t')] 9c ^ t y\ ( io) 

The final, nonlinear term of Eq . 1 1 Ul explicitly couples the 
dynamics of fluctuations at different wave vectors, so that 
the decay of C^ CA (t) is not simply exponential. The 
mode-coupling estimate of the static structure factor, 
C^ ICA (0) = l/p,™ CA , is determined by a self-consistent 
equation involving both the "tadpole" and "sunset" dia- 
grams. We avoid this static calculation by instead replac- 
ing /i^ CA in Eq. ^Jwith the exact form of C 1 [" 1 (0) from 
numerical simulations. This procedure is commonly em- 
ployed in mode-coupling studies of supercooled liquids. 

As the microphase transition point is approached from 
high temperature, we expect that only modes near the 



ordering wave vector ko will remain important. In this 
regime, a reduced model without reference to the cou- 
pling of specific length scales should capture the qualita- 
tive behavior of Eq. ^| Such a schematic model is sim- 
ilar to that studied by Leutheusser for structural glass- 
forming liquids [33 . (Indeed, Eq.^Jis only slightly differ- 
ent from that encountered in the idealized mode-coupling 
theory of supercooled liquids^]. In particular, the mem- 
ory kernel involves a two-point correlation function raised 
to the third, rather than second, power.) Restricting at- 
tention to |k| = ko, and neglecting coupling to other 
wave vectors, Eq. 1101 reduces to the dynamical equation 
exactly describing the p-spin model of a mean field spin 
glass (with p = 4). Since for p > 2 the critical properties 
of such models are essentially p- independent, we expect 
that near a critical point, Eo. 1101 will exhibit a plateau 
and an eventual transition to nonergodic behavior 38] . 

The character of the slow dynamics resulting from the 
theories underlying Eq. Q3] and Eq. ^| are fundamentally 
different. The dynamical Hartree theory (Eq.[8J) may ex- 
hibit a rapid slowing of dynamics as a function of inverse 
temperature only if the statics, as expressed through the 
renormalized mass th, are strongly temperature depen- 
dent. On the other hand, due to the nonlinearity of 
Eq. ^3 a slight change in the structure factor may re- 
sult in a dramatic change in relaxation times. It is well 
known that glass-forming liquids show little change in 
static structure as the glass transition is approached [13|. 
Thus, theories of the type given in Eq.[S]are not relevant 
near the glass transition. In the following sections, the 
predictions of Eq. [H] and Eq. ^| will be compared with 
simulations for the Coulomb-frustrated system. 

The renormalized perturbation theories developed in 
this section are strictly valid only in the limit of weak 
coupling, i.e., for small A or large, positive r. In our 
simulation work, we fix A = 1. It is thus instructive 
to ask at what value of r these theories are expected 
to break down. To answer this question, we follow the 
arguments of Hohenberg and Swift p4|. Specifically, we 
compare the Hartree approximation to the renormalized 
mass, th, with corrections introduced by mode coupling 
(i.e., the renormalized sunset diagram). These correc- 
tions are comparatively small when 



1-7-1 < 0.2 h 



7/5 



(11) 



In this regime the Hartree and mode coupling approxima- 
tions are controlled, and differ only quantitatively from 
one another. For larger |r|, however, the two approxima- 
tions can differ substantially, as we will see in numerical 
results presented in the following section. It has been 
noted previously that the static Hartree approximation 
can be accurate beyond its strict range of validity. There 
is thus no guarantee that a range of r exists in which 
the mode coupling approximation significantly improves 
upon an appropriately chosen harmonic reference system. 
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IV. SIMULATIONS 

In order to follow the dynamics of Eq. [2] or Eq. 0] nu- 
merically, it is necessary first to coarse-grain the field 
(f)(r, t) in space. This procedure yields a (periodically 
replicated) finite set of dynamical variables, whose time 
evolution may be integrated approximately over short in- 
tervals. We select a coarse-graining length a — 2ir/nko, 
and define new fields at lattice points r<: 

*i(t) =a- 3 f dr<j>(r,t), (12) 

where Vi is bounded by a cube of side length a centered at 
Yi. In the calculations described below, n — 8, so that a 
domain of wavelength ir/ko comprises several "soft spins" 
$i. To lowest order in a and a small time increment At, 
these renormalized fields evolve according to 

$i(t + At) = $;(i) - At[(r + fc - 2 (£ + fc ) 2 )$ 4 (t) 

+^m+m- (13) 

Here, the lattice approximation to the Laplacian op- 
erator acting on a function of space, Cf(xi) = 
a~ 2 ^2j(z nn f{xj) — f(%i), is taken to include a sum 
over nearest neighbors only (denoted nn). After coarse- 
graining, statistics of the random force remain Gaussian, 
with {rj(ri,t m )r;(rj,tn)) = 2(3' 1 (At/a 3 )S tj S mn . 

The simulation algorithm described above (which is 
very similar to those of Refs.[25l l3llp has several advan- 
tages over the numerical approach of Ref. |l8j , which em- 
ploys "hard spins" (<E>i = ±1) and explicit frustration. 
First, spins interact only with nearest and next-nearest 
neighbors, providing linear scaling of computational ef- 
fort with system size. Because cumbersome techniques 
associated with long-ranged forces are not required, a 
larger set of dynamical variables may be considered. In 
our calculations, the periodically replicated unit cell in- 
cludes 64 3 spins arranged on a cubic lattice. More impor- 
tantly, the coarse-graining procedure allows the dimen- 
sions of the unit cell to scale with the physically relevant 
length fc^ 1 . As a result, the unit cell spans several cor- 
relation lengths, even for very small values of fco- By 
contrast, in the work of Ref.jlg, t ne unit cell is com- 
parable to a single natural lamellar spacing for several 
of the simulated states (particularly those corresponding 
to "fragile" systems). In those cases, significant finite 
size effects are possible. In effect, Grousson et al. cut 
off slowly relaxing fluctuations at small k rather than 
the rapidly relaxing fluctuations at large k that are in- 
tegrated out in our approach. In Sec. 5 we discuss the 
dynamical implications of such a cutoff. 

Using Eq. 1131 we have computed the dynamics of sev- 
eral states of the model system at temperatures above 
the microphase separation transition (r > T tr ). We focus 
on three values of fco (0.1, 0.5 and 1.0) corresponding to 
a somewhat broader range of model parameters than was 
considered in Ref. [l8| . In each case, A = 1, so that the 




r < r tr (d). The configuration depicted in (c) was obtained 
from the nonequilibrium evolution of a nearly random state 
quenched instantaneously to r < rtr. 



microscopic dynamics is in principle strongly nonlinear. 
Representative configurations of the system are depicted 
in Fig. |21 typifying the high-temperature paramagnetic 
phase (a), the disordered phase near the microphase sep- 
aration transition (b), the ordered lamellar phase (c), and 
a nonequilibrium state produced by rapid quenching of a 
disordered system to low temperature (d). 

For each value of fco we consider, the relaxation of spin 
correlations slows dramatically near the transition to mi- 
crophase separation. The time required for single-spin 
correlation, 

f r/k 

C(t) ee ($i(0)Mt)) = N- 1 J ^C k (t), (14) 

to decay to 10% of its initial value, t, is plotted as a 
function of r in Fig. [3] The growth of relaxation times 
as r approaches r tr is sharpest for the smallest value of fco. 
Indeed, critical fluctuations are suppressed least strongly 
in this evidenced by the onset of sluggishness 

very near r = 0. 

Grousson et al. have likened systems corresponding 
to large and small values offcn to "strong" and "fragile" 
glass formers, respectively [l8|. For fragile cases, they 
have even shown that the temperature dependence of 
t is well fit by the Vogel-Fulcher form found for super- 
cooled liquids. Although this functional form suggests an 
emerging importance of activated processes, the dramatic 
growth of relaxation times is in fact well captured by the 
harmonic reference system described in Sec. 3. Numer- 
ical solutions of Eq. (plotted as solid lines in Fig. |3Jl , 
corresponding to this Hartree approximation, are espe- 
cially accurate in the most fragile case (fco = 0.1). Even 
for the least fragile case (fco = 1), computed rates dif- 
fer from predicted values by at most a factor of two. 



n I I i I i I i I i I i I i I i I I 

-0.3 -0.2 -0.1 0.1 0.2 0.3 0.4 
X 

FIG. 3: Time t required for spin correlations to decay to 10% 
of their initial values, as a function of scaled temperature 
t. Circles, squares, and diamonds show numerical results for 
fco = 0.1, fen = 0.5, and fco = 1.0, respectively. Solid lines 
are predictions of the Hartree approximation described in the 
text. Arrows indicate temperatures at which this approxi- 
mation is expected to break down for each value of fco, as 
estimated using Eq. 1111 



Activated barrier crossing is manifestly absent on a har- 
monic landscape, strongly implying that slow dynamics 
are driven by static renormalization, rather than by fun- 
damental changes in the structure of trajectory space[{|. 
The static structure factor, Ck(0), in fact becomes more 
sharply peaked in a way that mirrors the sudden growth 
in t. In other words, the slowing of relaxation appears to 
be an example of de Gennes narrowing |34j. 

The time dependence of spin correlations provides fur- 
ther evidence for this interpretation. Specifically, even 
when relaxation is very slow, correlations decay nearly 
exponentially at each wave vector. In Fig. Ck(t) is 
plotted for many k values for a system very near mi- 
crophase separation (fco = 0.5, r = —1). Included wave 
vectors span a range from the lowest accessible spatial 
frequency (fc = 2it / L) to several multiples of fco- In 
no case is relaxation detectably caged or stretched at 
long times. The single-spin correlation function, C(t), in 
Ecj-Elis a superposition of all Ck(i), and thus does not 
decay as a single exponential. At long times, however, re- 
laxation is dominated by the slowest modes, those with 
wave vectors lying in a spherical shell with |k| « fco, and 
is very nearly exponential. Since static correlations are 
strongest for these modes, especially near the microphase 
separation transition, C(t) is nonexponential only over a 
small range of the total decay. 

Remarkably, spin relaxation is essentially identical for 
a very different choice of microscopic propagation rules 
which conserve the field $. In this Monte Carlo dynam- 
ics, described in Sec. 2, each trial move simultaneously 
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FIG. 4: Ck(t) and In [Ck(i)/Ck(0)] for several |k| as a function 
of lnt and t, respectively, for the state fco = 0.5, r = —0.1. 
The uppermost curve in each panel corresponds to |k| = fco 
and exhibits the slowest decay. 

displaces $ at a randomly chosen site i and at a a site 
j randomly chosen from the nearest neighbors of i. The 
displacement at i, A$j, is exactly compensated by that 
at j, i.e., A$j = — A$j. In this way, $ is conserved at all 
length scales greater than or equal to the lattice spacing 
a. In general, such a constraint can influence dynamical 
behavior dramatically. For instance, scaling exponents 
for unstable domain growth in similar models depend 
intimately on the conservation of order parameters|33j. 
The relaxation described above, however, is modified by 
the constraint only at very short times. 

The insensitivity of slow dynamics to field conserva- 
tion in this model was anticipated by Sachdev, who noted 
that the corresponding constraint couples strongly only 
to fluctuations with very small wave vector (|k| « 01|10|. 
The long-lived correlations in Fig. however, are gov- 
erned not by these modes, but instead by fluctuations 
of finite wave vector (|k| « fco), which couple relatively 
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FIG. 5: Decay of spin correlation C(t) predicted by Hartree 
approximation (lines), compared with results of numerical 
simulations (symbols) for several thermodynamic states: r = 
-0.14, t = -0.12, r = -0.1, r = -0.08, r = -0.06, 
r = —0.04, r = —0.02, r — 0, (in order from top to bot- 
tom in the plot). In each case, ko — 0.5. 

weakly to the constraint. Conservation is thus of only 
modest importance at long times in the disordered phase. 
These facts can lead in principle to dramatic finite size 
effects in numerical simulations. Specifically, the smallest 
periodically replicated unit must accommodate fluctua- 
tions with wavelengths several times ko. Otherwise, the 
absence of truly long wavelength fluctuations will gener- 
ate spuriously strong coupling of conservation constraints 
to the slowest accessible modes. As a result, artificial dy- 
namical features may appear at long times. For model 
energetics that are similar to but different from Eq. ^ 
stretched exponential relaxation has been computed from 
simulations in which system dimensions are comparable 
to fee -1 03. To the extent that fluctuations near |k| = fco 
are independent of model details, such anomalous dy- 
namics should not survive in the thermodynamic limit. 

The persistence of exponential relaxation into the 
neighborhood of microphase separation, obtained numer- 
ically both for Langevin and Monte Carlo dynamics, is 
consistent with the dynamical Hartree approximation de- 
scribed in Sec. 3. In fact, this simple theory also predicts 
with remarkable accuracy the time scales of relaxation, 
even as they grow by several orders of magnitude. Re- 
sults for the case fco = 0.5 are compared in Fig. OD by 
plotting C(t)/C(0) obtained from theory and simulation 
for several values of r. 

The success of the harmonic reference system pro- 
vides further evidence that "vitrification" in this model 
is driven by dramatic changes in structural order rather 
than novel relaxation mechanisms. Indeed, several dy- 
namical features that would seem to be related to ac- 
tivated barrier crossing may be well rationalized in the 
Hartree picture. For example, the fragility parameter 



defined empirically by 

t oc exp ^ ^y^ (15) 

appears to scale as D ~ Q 1 ^ 2 in the simulations of 
Ref.|18|. Here, Tj< is a fitted Kauzmann temperature 
at which relaxation times appear to diverge. Schmalian 
and Wolynes have suggested that this scaling arises from 
the entropy of activation for structural rearrangement of 
mesoscopic domains 20, |2l|, |22j. The Hartree approxi- 
mation offers a simpler explanation. Fitted relaxation 
times follow t ~ r^ 1 ~ exp(D/r) rather well. It is easy 
to show by dimensional scaling that the renormalized 
mass is a function only of t/Q 1 / 3 , yielding immediately 
D ~ Q 1 / 3 !^]. This slightly different scaling form fits the 
results of Grousson et al. equally well 13 • For the simu- 
lations of larger systems we have presented, the Hartree 
prediction appears to be superior. 

The dynamical scenario predicted by the simplified 
mode coupling theory of Sec. 3, on the other hand, is 
not borne out in our simulations. Most significantly, we 
observe neither loss of ergodicity nor two-step relaxation 
over the relevant range of r > r tr . These failures of 
the dynamically nonlinear approximation are evident in 
Fig. H3 comparing the MCA and simulation results for 
the same system and thermodynamic states considered 
in Fig. 03 Interestingly, the MCA predicts trapping at 
values of r for which the Hartree approximation remains 
reasonable. The infinite series of terms incorporated in 
the mode coupling approximation thus adds little real- 
ism to the lower order description, and eventually leads 
to incorrectly anomalous behavior. This series of terms 
must be compensated to a large degree by omitted terms 
at each order. There have been suggestions that simi- 
lar cancellation occurs in mode coupling expansions of 
supercooled liquid dynamics js^]- In that case, however, 
signatures of idealized mode coupling (i.e., two-step re- 
laxation) survive despite the existence of omitted relax- 
ation channels. In our model, no such signatures are 
evident. 

We emphasize that the thermodynamic states we have 
simulated lie exclusively in the disordered phase of our 
model system. We have confirmed this fact by comput- 
ing the work to reversibly impose long-range order (i.e., 
nonzero (|3>k|) at |k| = fen). Although the static Hartree 
approximation suggests that the states of lowest tem- 
perature considered for each ko have global free energy 
minima in ordered configurations, the computed free en- 
ergy of the paramagnetic state is in fact lower for each 
case. This quantitative failure of Brazovskii's approxi- 
mation is not surprising, as the relevant states lie out- 
side the strict range of validity of the approximation for 
A = 1 [23, 01 • Ou r analysis thus leaves unexplored a nar- 
row window of the disordered state extremely close to 
t = T tr . While it is possible that qualitatively new dy- 
namical behavior arises in this region, it cannot account 
for the dramatic slowing down we have demonstrated at 
higher temperature, which is driven by extreme struc- 
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FIG. 6: Comparison of MCA (solid lines) with numerically 
simulated dynamics (symbols) (a) and MCA (solid lines) with 
Hartree approximation (dashed lines) (b). Results are shown 
for thermodynamic states identical to those plotted in Fig. 5. 
According to Eq. 1111 contributions of mode coupling become 
significant around r ~ —0.06 (fifth line from the top). Note 
that this temperature is very near the critical temperature at 
which MCA predicts loss of ergodicity. 



tural changes rather than activated processes. Further- 
more, the sudden onset of nontrivial behavior would be 
in stark contrast to the more gradual onset of caging and 
stretched exponential behavior in real liquids. 



V. CONCLUSIONS 

We have examined in detail a simple model that is 
closely related to the frustration-limited domain theory of 
Kivelson, Tarjus and coworkers^, and to the uniformly 
frustrated "stri p e- g lass" model studied by Schmalian 
and Wolvnes|2(l l2ll |22j. The disordered phase of this 
model system indeed displays some hallmarks of molec- 



ular glass-forming liquids. But its dynamics differ quali- 
tatively from generic glassy behavior in several respects. 



While relaxation times increase dramatically in these 
models in a non-Arrhenius fashion as temperature is low- 
ered, we find that an optimized harmonic reference sys- 
tem captures the time dependence of fluctuations semi- 
quantitatively. In contrast to the vitrification of molecu- 
lar liquids, the onset of this sluggishness is not accompa- 
nied by significant power law or stretched exponential de- 
cay of correlations in time. Perhaps most importantly, we 
find that the the slow decay of dynamical correlations is 
driven by significant changes in static structure. An ide- 
alized mode-coupling theory captures these changes less 
accurately than the simpler harmonic approach, predict- 
ing caging and eventual trapping at temperatures where 
the simulated dynamics remain exponential and ergodic. 
Together, these results strongly suggest that slow dynam- 
ics of this model system arise from the same Gaussian 
fluctuations that drive the microphase separation tran- 
sition first order. This conclusion differs markedly from 
those of Kivelson et a^QI. Grousson et a/.fisl. flS]. and 
of Schmalian and Wolvnes |2(l l2lL l22j . which invoke acti- 
vated barrier crossing to explain the dramatic tempera- 
ture dependence of relaxation times. Uniform frustration 
alone thus appears insufficient to account for the unusual 
relaxation properties of supercooled liquids. It remains 
possible, of course, that such frustration plays a key role, 
but in concert with other important physical mechanisms 
and constraints. 



The simulations described in this paper constitute 
a thorough study of the thermodynamically disordered 
phase of a uniformly frustrated model. They raise some 
interesting questions about such models and leave others 
unaddressed. In particular, while glassy behavior is not 
manifest for r ^ T tr , we can not rule out the existence 
of caging or trapping extremely close to or below the 
microphase separation temperature (since relaxation at 
these temperatures is prohibitively slow) . To explore this 
possibility, we have simulated several initially disordered 
states with r < r tr . Their nonequilibrium evolution to- 
wards long-ranged order exhibits the self-similar aging 
characteristic of coarsening phenomena, but for numeri- 
cally accessible time scales does not resemble the relax- 
ation of a supercooled liquid. Nonexponcntial equilib- 
rium relaxation could yet emerge at temperatures above 
but very near the order-disorder transition 39] . Eastwood 
and Wolynes have in fact recently suggested that, for cer- 
tain spin models, surface tension effects drive the onset 
of activated dynamics nearly to the ideal glass transition 
temperature^. Although quite interesting, such behav- 
ior would itself be uncharacteristic of glass-forming liq- 
uids, as it would appear only within an extremely narrow 
thermodynamic range. 
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